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Abstract 

The use of hyperasymptotics and the Weniger transformation has been proposed, in a joint 
fashion, for decoding the divergent asymptotic series generated by the steepest descent on a wide 
class of saddle-point integrals evaluated across Stokes sets [R. Borghi, Phys. Rev. E 78, 026703 
(2008)]. In the present sequel, the full development of the H-WT up to the second order in H is 
derived. Numerical experiments, carried out on several classes of saddle-point integrals, including 
the swallowtail diffraction catastrophe, show the effectiveness of the 2nd-level H-WT, in particular 
when the integrals are evaluated beyond the asymptotic realm. 
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I. INTRODUCTION 

The present paper is a sequel of a previous work[l| concerning the evaluation of saddle- 
point integrals of the form 

Iik) = j^gis)exp[-kf{s)]ds, (1) 

where C is a suitable integration path in the complex s-plane, g{s) and f{s) are functions 
which, for simplicity, will be assumed to be nonsingular, and k will be intended as a "large" 
(in modulus) complex parameter. As is well known, the numerical evaluation of integrals of 
the kind in Eq. ([1]) is customarily required for solving several classes of physical problems, 
occurring in optics, quantum mechanics, statistical physics, fluid mechanics, and so on. In 
optics, the evaluation of several diffraction integrals is customarily carried out asymptot- 
ically by identifying the parameter k as the wavenumber of the radiation|2|. In quantum 
mechanics, the same role is played by the inverse of the Planck's constant, while in fluid 
mechanics by the Reynold's number 3|. 

In the stationary phase treatment of diffraction integrals the values of the associated 
complex wavefield are asymptotically evaluated by taking the contributions coming from the 
stationary points of f{s), each of them associated to a "ray" in the corresponding geometrical 
picture. Of particular importance is the birth and the death, as the spatial parameters ruling 
the "phase" function i/(s) vary, of "evanescent" rays across sets of codimension 1, named 
"Stokes sets" Q,^. 

The 6-, or Weniger, transformation 0, Q, Q] (WT for short), is particularly efficient for 
resumming the factorial divergent asymptotic series well away from Stokes sets, as well as 
sets where two or .no.e saddles are synnnetncaUy p.aeed in the eomp.ex s.ngu.a,. spaeefl. 
Unfortunately, as with other resummation techniques [61. 110| . the WT fails to perform across 
Stokes sets. The reason of such a failure stems to the extreme " specialization" of the trans- 
formation itself, which requires, for a successful resummation, an alternating sign pattern of 



the sequence of the single terms of the series 111. Several methods have been conceived for 
resumming nonalternating, slowly convergent or divergent, sequences [8|, |12| . some of them 
being based on the serial combination of various resummation techniques jll 
class of saddle-point integrals in Eq. ([T]), the marriage between hyperasymptotics|l4j. Il5l| (H 
for short) and the WT generating the so-called H-WT (which stands for hyperasymptotic 
Weniger transformation), allows the WT to successfully operate also across Stokes sets. 
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Basically, the H-WT consists in the sequential application, to the integral in Eq. ([T]), of a 
classical hyperasymptotic treatment, as described in Ref. |l5|], followed by the action of the 
WT on all asymptotic divergent series generated by H. In particular, the results obtained 
have shown how the Ist-order H-WT, for which only the first-stage of H is anticipated to 
the WT, is able to provide relative errors several orders of magnitude smaller than those 
achievable via the use of full hyperasymptotic treatments and with considerably lighter 
computational complexity and effort. A key aspect is that, differently from H, the first 
truncation operated on the starting asymptotic series has not to be an optimal, in the sense 
of superasymptotics (i.e., at the least term) one, but rather the corresponding truncation 
order, say A^, must be used as a free parameter for the subsequent application of the WT. 

A question which was not addressed in Ref. but mentioned only in the last sentence, 
is whether WT and H can be combined to higher orders in H, and if so, how the accuracy 
improves with order. The present paper is aimed at giving a first answer to such a question. 
We shall limit our analysis only to the second stage of H. Further increasing of the H-WT 
order would be achievable along the same guidelines outlined here. On the other hand, 
it should be noted how, on increasing the order of H, the number of asymptotic series 
associated to the corresponding remainder that have to be resummed exponentially grows 
for topologies involving more than two saddles and, at the same time, the number of free 
parameters (i.e., the truncation orders at each H step) linearly increases. Accordingly, from 
a mere computational viewpoint, it is mandatory to find a compromise between the H-WT 
order and the computational effort. Some limits of the Ist-order H-WT have already been 
emphasized in Ref. ]|, where asymptotic evaluations of saddle-point integrals for "small" 
values of the asymptotic parameter were considered. In such cases, in fact, to increase 
the parameter A^ of the Ist-order H-WT does not necessarily lead to an improvement of 
the reached accuracy, but often results in the opposite, i.e., a worsening of it. It is just 
the above scenario that we are interested in when the 2nd-level H-WT will be developed. 
Numerical experiments will be carried out on the class of saddle-point integrals already 
considered in the numerical sections of the first paper |l| . Moreover, asymptotic evaluations 
of the so-called swallowtail diffraction catastrophe will be proposed as a new numerical 
experiment. The swallowtail function is defined via Eq. ([1]) with f{s) being a 5th-order 
polynomial with respect to s, thus involving a four saddle network. We will present a study 
of the accuracy achievable via H-WT asymptotic evaluations of the swallowtail diffraction 
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catastrophe for points placed at the Stokes set, following the prescriptions by Berry and 
Howlsjsl. In doing so, we will find that the corresponding asymptotic expanding coefficients 
can be expressed in closed-form terms. 

For any practical implementation of the H-WT, a key role is played by the numerical 
evaluation of the corresponding hyperterminants[l4. [l5| which are defined through suitable 
multiple integrals. For the lowest-order hyperterminant the exact analytical expression is 
available from literature [l^, but unfortunately this is not true for higher-order hyperter- 
minants, included those involved in the 2nd-level H-WT. In the present paper we solve the 
problem of the 2nd-level hyperterminant exact evaluation for a particular, but very impor- 
tant, choice of the hyperterminant parameters, which often occurs in the implementation of 
H for evaluating a wide class of saddle-point integrals. Up to our knowledge, this is a new 
result which also provides an interesting connection of such hypeterminants to the Meijer-G 
functions 17|]. Moreover, although the closed- form evaluation of 2nd- level hypeterminants 



for arbitrary choices of their parameters seems to remain an open problem, in the present 
paper we find a semi-analytical representation which turns out to be suitable for numerical 
calculations via standard integration packages. Similarly as done in Ref. jll, one of our 
aims is to keep the paper reasonably self-consistent. Accordingly, in the next section a brief 
review of H, up to the 2nd-level, is given. As far as the WT is concerned, we believe that 
what is contained in Ref. jl], together with the extensive bibliography, should be enough 
also for a nonexpert reader. For this reason, we do not repeat it in the present paper. 



II. RESUMING HYPERASYMPTOTICS 
A. Preliminaries and notations 



For simplicity, we shall refer to the asymptotic evaluation of saddle-point integrals of the 
type in Eq. ([1]) where the set of saddle points of /(s) will be denoted S and the integration 
path C will be thought of as the union of a finite number of steepest descent arcs each of 
them, say C„, passing through the contributive saddle point s„, which will be supposed to 
be a simple one. Accordingly, the quantity can generally be written as 

I{k) = I g{s) eM-kf{s)]<ls=Y, ^^"H^), (2) 
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where S' denotes the subset of S containing all the contributive saddles, and 

J(")(fc)= / gis)exp[-kfis)]ds. (3) 

The last integral can be written asfisl 

Ir^ik) = k-^l^ eM-kfn) T^'^Kk), (4) 

where /„ = f{sn), and where T^'^\k) can formally be written through the following asymp- 
totic series expansion: 

oo 

T(")(A;) = ^ A;-"TJ"), (5) 

r=0 

the expanding coefficients tJ"^ being expressed via the integral representation 

2vri I [/(.)- W 
where the subscript n denotes a small positive loop around the saddle s„. 

B. Development of H up to the second stage 

H starts by writing Eq. (|5]) in the form 

7V-1 

T(") {k) = Y^ k-'' T/") + {k,N), (7) 

r=0 

where N represents a positive integer and R''^\k,N) = XI^at k~^ Tr'^\ denotes the corre- 
sponding remainder which, due to the diverging character of the asymptotic series, turns 
out to be a diverging quantity too. H is based on a couple of fundamental results, found via 
a nontrivial analysis in Ref. 15(]. The first is that the value of the expanding coefficients 



Tr"^ at the saddle s„ is strictly related to the values of the expanding coefficients Tr"^^ at 
all those saddles, say {sm}, which are adjacent to s„, via the following formal resurgence 



relation 



Ti'" = ^ E (-1)"" E (8) 

where An denotes the set containing the indexes pertinent to all saddles adjacent to s„, the 
quantities Fnm, called singulants, are defined by 

Fnm fm fm (9) 



and the binary quantities 'jnm ^ {0, 1} are obtained through a topological rule[l5l|. The other 
fundamental tool of H is the following integral representation of the remainder i?^"^ [k, N) 

1 ^ (_l)7nm 



3: 



^ ' ^ 27ri ^ (kF, 



N 

nm j 



(10) 



kF 

Equations (|3])- (ITU]) allow hyperasymptotic expansions for the saddle integral in Eq. ([2]) to 
be built up in principle to any order [l^. For instance, the direct substitution of Eq. ([5]) 
into Eq. ffTOl) leads to the Ist-stage hyperasymptotic expansion 



r=0 



m<=An 



oo 

X 

T- = 



where the function Kn\p), called hyperterminant of order iQ, [3], is defined through the 
integral 



^Wrm 1 Ha exp(-t;) 

/5 



where, in order for it to converge, n > 0. Moreover, it can be shown that 



K^:\^) = exp(/5) ^ in - 1)! + (-1)"-M7re exp(/?), (13) 



where -En(') denotes the exponential integral function [17|, while e equals 1 if /? < and 
zero otherwise. The presence of the term containing e has to be ascribed to the evaluation 
of the integral in Eq. (IT^ . when /? < 0, in the Cauchy principal- value sense. Equation (fTTj) 
represents the first hyperasymptotic stage, at which the divergence of the asymptotic series 
in Eq. (jS]) is led back to the presence of adjacent saddles 1^. Furthermore, the asymptotic 
series in Eq. (fTTj) are only formal, since for r > N the terminant -ft'^Lr diverges. In Ref. l| 
Eq. (|TT1) was taken as the starting point for introducing the H-WT. In particular, instead 



of using the WT directly on the single terms of the series in Eq. ([5]), it is employed for 
resumming the asymptotic series associated to all saddles Sm, with m G An, which appear 
in Eq. fill I) . Of course, due to the fact that r < in Eq. (ITT]) , it is mandatory that be left 
as a free parameter, in order for the WT to be able in decoding the above asymptotic series. 
The 2nd-level H can be derived by truncating each of the asymptotic series in Eq. fill I) at 
an order, say M, and by generating, for each adjacent saddle Sm, with m G An, a list of 
asymptotic series associated to all saddles, say Sh, such that h G Am- In Appendix |Xl only 
for the reader convenience, the derivation of the 2nd-level hyperasymptotic expansion of the 
integral in Eq. ([T]) is briefly recalled according to the formalism of Ref . [l^ . In particular, 



it is found that 



N-l 



r=0 



2vri 



M-1 



X > {-lyk-'-Tl^^K^ZA-kFnm) (14) 



r=0 



(27ri)2 



E E (-1) 



7nm+7mh 



OO ^ 

V A'-^T^'^) /r^^) -kF ■ 



oo 

X 

r=0 



mh 



'(2) 



where Kn^m{(3; 7), the hyperterminant of order 2, is now defined through the double integral 



00 /"OO ^nn-l^n-l 



(15) 



X / / dti df exp(— u — 7t> 
'0 Jo 



Similarly as we done in Ref. pj], Eq. fJHj) can be used to give estimates of T^^\k), as 
functions of the two (free) parameters A^ and M, by resumming, via the WT, all asymptotic 
series generated at the second stage of H which are inside the double sum with respect to h 
and m. 
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III. ON THE EVALUATION OF 7) 

The numerical evaluation of the hyperterminants represents a step of fundamental impor- 
tance for any practical implementation of the H-WT algorithm. Unfortunately, differently 
from the lowest-order H-WT, for which the corresponding hyperterminants are achievable 
via the closed-form expression in Eq. f|T3|) . there are no analytical expressions available for 
higher-order hyperterminants. In a series of important papers, Olde Daalhuis 
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ad- 



dressed the general problem of the hyperterminants evaluation, up to arbitrary precisions, 
through the use of convergent series representations based on hypergeometric functions. 
However, for the particular case of the 2nd- level hyperterminant, it seems that some new, 
at least up to our knowledge, results could be established. 

From Eq. (fT5|) where, in order for it to converge, n > and m > 0, the hyperterminant 
can be written as 

,(2) . 1 



(16) 



foo i>co ^m^n-l 



poo poo 

X / d-udf exp(— M — 7t>) 
Jo Jo 



1 + ^1 (u + v) 



and, by formally expanding the factor 1/(1 as a geometric series, after some algebra 



takes the form 



^-(ft,)^<=iPi:(-i) 

k—m ^ ' 



00 /"OO yk^n-\ 



X / d-u df exp(— M — 7 ti) 
Jo Jo 



'-V. 



u + v 
n - 1)! 



(17) 



^ / 1\' fc! / , 1 
^ I] F ( n, 1; /c + n + 1; 1 - 



k=m, 



(3 J k + n \ 7^ 
where F(-, ■; ■; ■) denotes the hypergeometric function Although the series in Eq. ( IT7|l 



is divergent, it can be decoded via Borel summation, i.e., by replacing the term k\ by its 
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integral representation, i.e., 

I'OO 

k\= dt exp(-^)^^ 
Jo 

which, once substituted into Eq. f|T71) . leads to 

Xa(/3;7) = (-l)'"^^^ exp(-t) 



°^ ' tV 1 



(18) 



(19) 



k=m 



F (n,l-k + n + l-l--] . 

p J k + n \ 7^ 



Although, as we shall in a moment, it is possible to express the series inside the last equation 
through a closed form, it is better to carry out the evaluations for the case 7 = 1 and 7 7^ 1 
separately. 



On letting into Eq. f[T9|) 7 = 1 we have 



K^uUP; 1) = (-1)' 



(n- 1)! 



(20) 



(n- 1)! 

(m + n)l3'^-^ 



m + n, l;m + n + l;— — ) . 



The integral in Eq. ( |20l) can be evaluated by using the representation of the hypergeometric 
function given by formula 9.34.7 of Ref. 17|]. In particular, it turns out that 



(n-l] 



n—l 



x/^ dtexp(-t)Q 



m+l 



'-^22 



(21) 



—n — 171,-1 
— 1, — n — m — 1 



where G^I^{-) denotes the Meijer function [iTj]. Finally, by using formulas 9.31.5, 7.813.1, 



and 9.31.2 of p/7[, after some algebra it is found that 



1 — n — m, 1 

1 — 11,1—11 — 171,0 



(22) 



Equation (122|) represents one of the main results of the present paper. As we shall see in the 
numerical section, in applying the 2nd-level H-WT the evaluation of the hyperterminants 
Kn}n{(3;'y) is often required for 7 = 1. This happens whenever the contributive saddle s„ 
turns out to be adjacent to itself after two hyperasymptotic stages, i.e., when h = n into 
Eq. mn 

For 7 7^ 1, the series inside the integral in Eq. (ITTI) can still be expressed through a 
closed form, although the subsequent integral unfortunately not. However, a semi- analytical 
expression, which turns out to be suitable for being evaluated via standard numerical inte- 
gration packages can be derived. In doing this, Eq. (ITTI) is first rewritten as 

. (n-l)\ 



X 



k=0 



P J k + n 



k\ ( , 7-1 
F I n, 1; + n + 1; 

7 



(23) 



where 



■5 = 1: 



fc=0 



1 



j3 J k + n 



B f T — 1 

F [n,l;k + n + l;- — 

7 



(24) 



so that the task is to evaluate the series in Eq. fl2^ for 7 7^ 1. On substituting from Eq. ([18 
into Eq. fl2^ we have 



/>oo 

S = I dt exp(-t) 
Jo 



X 



E 

A:=0 



k F n, l;A; + n + l;l 

_i\ _v 7 

(3) k + n 



7 / dt exp(— t) 
Jo 



(25) 



X 



E 

fc=0 



t\ F(l + A;,l;n + l + fc;l-7) 



k + n 



where use has been made of the relationfsee Ref. 



20| p. 347] 



F n, 1: A; + n + 1: 1 



7 



7F(l + fc, l;n+l + A;;l-7). 



(26) 
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Finally, on writing Eq. (125|) as 



s = l 

n 



POO 

/ dt exp(— t) 
Jo 



X 



V (1 



(27) 



fc=0 



- F{l + k,l;n + l + k-l--f), 



where 



)k denotes the Pochhammer symbol, formula 6.7.1.8 of Ref. 20| gives at once 



s = l 

n 







1 + 



t 



\ 



l,l;n + l:l 



7 



V 



1 + 



t 



(28) 



Pj 



Notice that, when Re[/5] > 0, the function S can also be evaluated through the alternative 
form 



S 



/37 



n 



1 exp 
dp 



-P 



1 







P 



F(l,l;n+l;l-7p). 



(29) 



Although it seems that the above expressions cannot be further simplified, the numerical 
evaluation of the function S can be done with high accuracies by using standard integration 
packages. Finally, it should be stressed that, for (3 < 0, the evaluation of the double integral 

(2) 

defining Kn,m{P]'y) has to be done, with respect to the t>-variable, in the Cauchy principal 
value sense, in order to overcome the singularity placed at f = —f3. This, in turn, implies 
that an extra term must be added to the result. In Appendix [B] such term is analytically 
evaluated starting from the definition in Eq. (fT5|) . and turns out to be 



(30) 



All subsequent numerical experiments will be done within the Mathematica language. 



IV. NUMERICAL EXPERIMENTS 



A. Evaluation of the Airy function across the Stokes Une 



Consider the evaluation of the Airy function, defined as 



Ai(a;) = 7^ / exp 



1\ — +XS 



ds, 



(31) 



11 



which is of the form given in Eq. with g{s) = l/(27r), f{s) = — i(s^/3 + xs), and k = 1. 
The detailed analysis of the saddle topology, as well as the expanding coefficients tJ"'* has 
been summarized in Ref. so that here it will not be given again. We only recall that 
the two saddles are Si = (— x)^/^ and S2 = — Si, and that Ai = {^}, 712 = 0, 721 = 



We focus our attention on the evaluation of the Airy function across the Stokes line 2l|] , 



i.e., for arg{x} = 2it/3, in order to compare the performances of the 2nd-level H-WT with 
respect those displayed by the Ist-order H-WT in the same situation. More precisely, we 
write the argument of the Airy function (3/4 X exp(i27r/3), where F IS a 

real positive parameter, whose value coincides with the singulant F12. The study of the 
asymptotic evaluation of the Airy function across its Stokes line has played a pivotal role in 
the development of several asymptotic technique, mainly in light of the relative simplicity 
of the involved saddle topology. Such a simplicity could help in grasping, whereas possible, 
some conceptual aspects related to the use of the H-WT. Differently from what done in 
Ref. where the relative error values were displayed via the use of tables, in the present 
paper we are going to resort to graphical visualizations, due to the presence of the two "free" 
parameters and M. In the first experiment, whose results are shown in Figure [Tj the Airy 



function is evaluated for F = 16. Note that the same experiment was carried out in Ref. [14 1 
via the use of H. The values of the relative error, obtained through the Ist-level H-WT, are 
shown, as black dots, versus the values of A^, reported on the abscissa axis. For each value 
of A^, the values of the relative error obtained via the 2nd-level H-WT, with M G [3, A^ — 1], 
are also plotted and, for the sake of clarity, are joined with lines of different color each of 
them, which corresponds to a different value of A^, departing from A^ itself. This can be 
noted from the figure, where it is immediately seen how the higher the A^, the longer the 
corresponding coloured "leg" is. From a first look to the figure, it appears that the relative 
error, obtained with both the 1st- and the 2nd-level H-WT, is lower bounded. We shall find 
that all subsequent numerical experiments present the same characteristic. As a first remark, 
it should be noting that the improvement of the estimate accuracy induced by the 2nd-level 
H-WT with respect to that obtained via the Ist-level one of the same order appears, at 
least for values of A^ not too large, not to adequately refund the unavoidable increase of the 
computational complexity required by the application of the 2nd-level transformation. For 
example, it is seen from Fig. [T]that a relative error of the order of 10~^^, achieved through 
the 2nd-level H-WT with A^ = 11 and M = 8, would be reached via a Ist-order H-WT by 
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letting = 13, but with a considerable saving of computational effort. The above example 
clearly suggests the use of the 2nd-level H-WT only in those cases where the best accuracy 
attainable via the Ist-level H-WT turns out to be not adequate. This happens, for instance, 
when the integral is attempted to be evaluated beyond the asymptotic realm. To put into 
evidence this aspect. Figure [2] shows the same as in Fig. [H but for a decreasing sequence 
of values of F, namely 14 (a), 10 (b), 6 (c), and 2 (d). In particular, in Fig. Myd), where 
the Airy function argument is located at a distance ~ 1 from the origin of the complex 
plane, the Ist-level H-WT provides a best error of the order of 10~^, achieved for = 4. 
Higher accuracies are not allowed because the information gained at the first H stage turns 
out to be no longer sufficient to generate WT-resummable sequences. The 2nd-level H- 
WT, on the other hand, provides a best error of the order of 10~^°, which is attained for 
{N,M) = (15,11). Some intuitive insights about the resummation process associated to 
the H-WT could be grasped by noting that a lower bounded error is an intrinsic imprint of 
superasymptotic and hyperasymptotic resummations, whereas it is not generally featured 
by the application of the WT to alternating factorial divergent series js]. Accordingly, one 
should be inclined to think that such an error behavior could be ascribed to the presence of 
the "regularization" step operated by H on the raw input data. Speaking within a more 
general context, this should be somewhat related to the possible presence on nonanalytic, 
nonperturbative correction terms which cannot be grasped simply by resummation processes, 
but rather require the use of "generalized nonanalytic expansions" Q]. 
In a second experiment concerning the Airy function, the asymptotic parameter F is let 
running within the interval [2,4] and, for each value of F, an exhaustive search of the 
optimal values of the truncations and {N,M), which minimize the 1st- and 2nd- level 
relative errors, respectively, is operated. The results are shown in Fig. [3l where the optimal 
relative errors obtained via the 1st- (open circles) and the 2nd-level (dots) H-WT are shown 
as functions of F. The values of the optimal truncation for the Ist-level H-WT are also 
reported, versus F, in Fig. HI while those of A^ and M, for the 2nd-level H-WT, in Fig. [5]^a) 
and (b), respectively. We will come back later on the above results. 

In concluding the present section, however, we want to provide a table of explicit values 
obtained through the use of the 2nd-level H-WT. We choose to evaluate the Airy function 
for F = 2, for which the optimal setting of the truncation parameters turns out to be 
{N,M) = (15,11). The preliminary step is the evaluation, through a simple WT, of the 
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contribution, to the Airy integral, coming from the saddle S2- The result is shown in Table 
[H Furthermore, the subsequent action of the 2nd-level H-WT on the saddle si is shown in 
Table [TTl where the complete estimates of the Airy function provided are reported together 
with the corresponding values of the truncation M G [3, 14], with = 15. 



B. Instanton integral 

The second numerical experiment concerns the evaluation of the instanton integral 

/ + 00 
exp[-k{s^ - If] ds, (32) 
-oo 

with k > 0. al^eady co,.ide.ed in Ref. Q as a sin^pMed prototype fa the ,^deU„g of 

instanton tunneling between symmetric double wells. It was shown in Ref. [IJ that the 
integral in Eq. fl32|) can be written as[l| 

Afik) =2k^/^Re{I{k)}, (33) 

where, by referring to Eq. ([T]), g{s) = 1, /(s) = (s^ — 1)^, and where C is the steepest descent 
path connecting the points — ioo and +oo via the lines Im{s} < and Re{s} > 0. The 
complete saddle topology, as well as the expressions of the expanding coefficients associated 
to all saddles have been described in Ref. [l|. In particular, there are three saddles, si = —1, 
S2 = 0, and S3 = 1, with Ai = {2}, A2 = {1, 3}, and ^3 = {2}. The saddles involved in the 
evaluation of X(/c) are S2 and S3, but only the latter requires a H-WT treatment, since the 
associated singulant is F32 = 1 > 0, and the corresponding asymptotic series turns out to 
be nonalternating. Furthermore, 712 = 1, 721 = 723 = 0, and 732 = 1, while we recall that 
the integral in Eq. fl32|) can be expressed in closed form via 



^f(k) = exp(-A;/2) 



k\ . fk 



(34) 



where /«(■) denotes the rath-order modified Bessel function of the first kind. The first 
experiment concerns the evaluation of A/'(l/2) via the 1st- and the 2nd- level H-WT. In 
Fig. E] it is seen how the 2nd-level relative error is bounded, with a minimum value of the 
order of 10~^, achieved for (A^, M) = (6, 5). On the opposite, the Ist-level H-WT turns out 
to be completely inadequate to provide a reasonably accurate estimate of the function, due 
to the very low value of k. The searching for optimal values has also been carried out in 
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the present case, but using, as the varying asymptotic parameter, k G [1/2,3], i.e., where 
the Ist-order H-WT displays the worst results in terms of accuracy, as shown in Fig. 2a of 
Ref. [l]. The error values are shown in Fig. [TJ versus k, while the optimal settings of and 
of (A^, M) are plotted, against k, in Fig. [S] and [HI respectively. 

It is worth now comparing the results pertinent to the Airy and the Af{k) functions. 
What we are going to show can seem at first sight somewhat surprising, but gives a possible 
first hint toward the understanding of the H-WT mechanisms. For simplicity we shall refer 
to the Ist-level transformation, but the results will apply also to higher-order levels. In 
Fig. [ini the values of the relative error obtained for the Airy function (dots) and for the 
instanton function (solid curve) are plotted, versus A^, when the values of the parameter 
F and the parameter k are numerically equal. In particular, figure (a) corresponds to 
k = F = 3, {h) to 7, (c) to 12, and (d) to 20. It is clearly seen that the behavior of the 
relative error follows basically the same law. To give a possible explanation of this, in Fig. [11] 
a pictorial representation of the complete saddle network and the complex integration path 
involved in the evaluation of the Airy (a) and of the instanton (b) functions is plotted. In 
both pictures, the black dot denotes the saddle for which the H-WT is required. Although 
the two saddle distributions are clearly different, they present some common features that, 
together with Eq. (fTTj) . are enough to justify what happens in Fig. [TOl Each of the "black" 
saddles is adjacent to a single saddle. For the Airy function si is adjacent to S2, while for 
the instanton function S3 is adjacent to S2- The values of the corresponding singulants are 
F and 1, respectively. The use of the resurgence relation in Eq. (IHD now gives, for the two 
"black" saddles, 

Ti" cc (35) 

for the Airy function and 

« (36) 

for the instanton function. From the above equation it is seen that the behavior of the 
expanding coefficients follows the same asymptotic law as soon as F = k. At the same time, 
however, the above equality guarantees that also the asymptotic laws for the expanding 
coefficients corresponding to the adjacent saddles is identical. In fact, for the Airy function 
the saddle adjacent to S2 is Si itself, with a singulant value of —F. As far as the instanton 
function is concerned, the saddles adjacent to S2 are Si and S3, but for both of them the 
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singulant values are -1. Accordingly, the use of Eq. ([8]), together with the condition F = k, 

(2) 

provides again an equivalence between the asymptotic laws of for the Airy function and 
Tr^'* for the instanton function. Finally, on using Eq. fill I) it is not difficult to convince 
that the retrieving process is the same for the two functions at the lst-level[23]. Leaving a 
deeper understanding of this phenomenon to future investigations, it is here worthwhile to 
point out how an immediate consequence of the above described "topological equivalence" 
could be the restriction of the study of the H-WT retrieving performances to a few classes 
of prototype test cases. 



C. Swallowtail diffraction catastrophe 



As a last numerical experiment we consider asymptotic evaluations of the so-called swal 
lowtail diffraction catastrophe 
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24 



25 



S{x,y,z) 



exp 



26l | , which is defined through the following integral: 

ds, (37) 






X ■ 



y^ + zs 



5 3 

which is of the form given in Eq. ([T]) with g{s) = 1, f{s) = — i(s^/5 + xs^/3 + ys'^/2 + zs), 
and k = 1. The integration path C can be thought of as the union of steepest descent paths 
approaching, for |s| ^ 1, the directions = {2n + l/2)7r/5, with n = 0, 1, ...,4. Although a 
systematic treatment of the swallowtail asymptotics, along the general classical rules recalled 



in Sec. Ill Al can be derived by paralleling the analysis carried out, for the Pearcey function, 
in Ref. [l5|, up to our knowledge it is not present in the current literature. As shown in 
appendix [Cl all expanding coefficients tJ""* are given by 

(5i)r.+l/2 1/2)1 



rpin) 



where 



(10s3 + 5SnX + 5?//2)5''/3+l/ 



2 ^2r ' 



a, (3), 



(38) 



a 



and where the polynomials B^\u,v) are defined via the generating function formula 

oo ^ 
n=0 



(104 + 5s„a; + 5i//2)V3^ 

lOs^ + 5x/3 
(10s3 +5s„x + 5?//2)2/3' 



(39) 



{t^ + ut^ + vt + l)^' 



(40) 
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It is not difficult to prove that Eq. (HO!) allows the numerical evaluation of the polynomials 
Bi^\u,v) to be efficiently performed via the use of the following recurrence rule, whose 
derivation is outlined in Appendix [D] 



nB, 



(A) 



n-3 + 3X) Bi% -u{n-2 + 2X) Bi% 



-vin~l + X)Bi'\, 



>(A), 



1,b['\u,v) 



-Af , and B. 



(A), 



U,V} 



(41) 



with the triggering values Bq 
l)/2. 

The numerical experiments we are going to illustrate concern asymptotic evaluations of 
S{x,y,z) at points belonging to the corresponding Stokes set, which has been estensively 



studied in Ref. 



51] (see, in particular. Fig. 3 of this reference). Accordingly, the triplets 



5| , in order to investigate 



(x, y, z) have been chosen following the prescriptions given in Ref. 
points at the intersection between the Stokes surface and the plane x = 0, along the branch 
corresponding to ?/ > and 2; > 0. This leads to triplets of the form (x, z) = (0, k^/^, x 
0.23012 . . .), with k being a positive parameter Q]. We start by considering the case k = 2. 
The saddle topology is constituted by four saddles, which are listed below together with the 
corresponding list of adjacent ones: 

si = 0.8062... -i 1.2357..., ^1 = {4}, 



S2 



S3 



'15 



-1.2828 



A2 = {4}, 



^3 = {4}, 



(42) 



S4 = -0.3296. 
with the orientation matrix being 



A = {1,2,3}, 



(43) 



1 1 ■ 

where the values of 'jnm which are not relevant for the present experiment have been replaced 
by dots. Three of the four saddles, namely S2, S3 and S4, do contribute to the integral. In 
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particular, the integration path consists in the union of three steepest-descent arcs, the 
first connecting the points oo exp(i97r/10) and oo exp(il37r/10) passing through S3, the 
second connecting the point 00 exp(il37r/10) to the saddle S2, passing through S4, and 
the third connecting the saddle S2 to the point 00 exp(i7r/10). Accordingly, the Stokes 
phenomenon occurs via the so-called saddle connection between saddles S4 and S2^- A 
pictorial representation of the topology above described is given, for the reader convenience, 
in Fig. [121 

Before showing the numerical results about the performances of the 1st- and the 2nd-level 
H-WT, it is worth giving some details about the way the saddle topology of the swallowtail 
integral influences the retrieving capabilities of the WT. We refer, in particular, to the 
contribution, to the swallowtail integral, of S3 and S4. As far as the former is concerned, 
since there is only one adjacent saddle, S4, due to the fact that the corresponding singulant 
-F34 = i 0.602168... is purely imaginary, it turns out that Tr^^ oc (— i)*" (r — 1)!/|F3^4|'', 
thus allowing the WT to operate the resummation. Similar considerations can be done for 
the contributive saddle S2. The situation is somewhat different for the saddle S4, which is 
connected to S2. In fact, from the above described topology, it turns out that the adjacent 
saddle S3 is dominant, i.e., presents the minimum value of |-F4,m|, with m G ^4. Accordingly, 
one should conclude also for Tr^"* an asymptotic "factorial divided by power" law, similar 
to that corresponding to Tr and, due to the fact that F4 3 = —i 0.602168. . ., one should 
expect the WT to be able in resumming the corresponding asymptotic series. But this, on the 
contrary, does not happen, because of the presence of the other two, nondominant, saddles 
si and S2, symmetrically placed in the complex singulant space. This can be explained by 
expliciting Eq. ([8]) as 



Stokes sets do not necessarily display a strictly nonalternating sign pattern as, for example, 
happened for the Airy function, but rather how the asymptotic behavior of their single terms 
can display more complex patterns, depending on the whole saddle topology. Figure [T3] 
shows the relative errors obtained through the 1st- and the 2nd-level H-WTs in the case of 
swallowtail evaluation across the Stokes set defined above, when k = 2. As for Figs. [H [2l 
and ini the errors are plotted versus the parameters N and M. In this case their optimal 




(44) 



where Tl 



{[Tq ]*. This example shows how the divergent asymptotic series generated on 
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values turn out to be = 4 (for the Ist-level) and (A^, M) = (7, 6) (for the 2nd- level), with 
corresponding error values of 2 ■ 10~^ and 5 ■ 10~^, respectively, evaluated with respect to the 
"exact" value of S{0, 2.8284 . . . , 0.9205 . . .), obtained via the method recently proposed in 
Ref. t\. Finally, an experiment about optimal resummation of swallowtail functions has 
been carried out, by using k G [2,4] as the parameter representative of the "asymptoticity" 
features. The error values are shown, versus k, in Fig. [TU while the optimal settings of 
and of (A^, M) are plotted, against k, in Fig. [15] and [161 respectively. 

V. CONCLUSIONS 

The H-WT was introduced in Ref. l| as a powerful and easily implementable refinement 
of H aimed at allowing the WT to successfully decode those divergent asymptotic series 
generated through the application of the steepest descent method to saddle point integrals 
evaluated across Stokes sets, for which their single terms do not display a strictly alternating 



sign pattern. The scheme proposed in |l| employed the WT only on the asymptotic series 
generated by the first-stage hyperasymptotic treatment of the corresponding diverging re- 
mainder. In the present sequel we reported on the possibility of combining H and the WT 
to higher orders in H. In particular, the full development of the 2nd-level H-WT has been 
detailed within the classical framework of H. The results obtained from the application of 
the 2nd-level H-WT, also in comparison to those obtained via the Ist-level one, on the dif- 
ferent types of saddle-point integrals considered, showed how the increase of complexity and 
computational effort required by the new transformation be adequately refunded, in terms 
of accuracy of the estimate, particularly when the integrals are evaluated for values of their 
parameters which are beyond the asymptotic regime, whereas H turns out to be inapplicable 
and the Ist-order H-WT unavoidably lacks precision. At the same time, however, it should 
be noted how, for "ordinary" asymptotic evaluations, at least in the cases considered in the 
present work, the performances of the 1st- and the 2nd-level H-WTs seem to be comparable 
in terms of the estimate accuracy, against a considerably difference in the computational 
efforts required by the two transformations. 

Although the H-WT has been developed, here and in Ref. []]], with reference to the 
evaluation of saddle-point integrals of the type in Eq. ([T]), we believe it could be useful also 
in dealing with problems of different nature like, for instance, the hyperasymptotic treatment 
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of a wide class of linear and nonlinear ordinary differential equations, which has been recently 



considered [28|, |29|, |30| . The semi-analytical algorithms proposed for the numerical evaluation 
of the 2nd-level hyperterminants would reveal useful in this perspective. 

Before concluding, it is worth pointing out, as an important open problem, the need for an 
a priori algorithm for estimating the values of and M that lead to optimal results. Espe- 
cially in cases where it is not convenient (or possible) to evaluate the original function, such 
an algorithm would certainly be of great help to a typical user. Unfortunately, differently 
from H, for which the optimal settings of the hyperseries truncations are directly extracted 
from the singulant values 0, Q, at present it does not seem possible to provide similar 
information for the H-WT. The difficulties in giving practical guidelines for the choice of 
and {N, M) can also be appreciated from the results presented in Sec. [IV] and especially 
from Figs. HI O [H [H [15], and [TH where it seems quite difficult to obtain general rules for the 
optimal settings of them. A possible hint, grasped from a quantitative comparison between 
the results obtained for the Airy and the instanton functions, seems to be given by the strong 
connection between the H-WT retrieving performances and the saddle topology associated 
to the integral under consideration. What we found is that different saddle networks can 
share a sort of "topological equivalence" property, which is related to the set of the saddles 
adjacent to that under consideration and to the values of the relevant singulants. If two 
networks turn out to be equivalent at a certain hyperasymptotic level, this would result in 
the same computational effort, in terms of relative error, as far as the corresponding H-WT 
retrieved estimates are concerned. This, in particular, would imply that the study of the 
H-WT retrieving performances could be, in principle, carried out only on a restricted class 
of prototype functions. 
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APPENDIX A: SECOND-STAGE HYPERASYMPTOTICS 

The starting point is again the integral representation of the remainder R^^\k, N) 
in Eq. flTUl) . where we let 

A/-1 oo 



(m) 

with M being a new truncation order. Substitution of Eq. (lAip into Eq. (fTOj) gives 



9'n-i ' ^ 



2vri 



A/-1 



r=0 



+— V (-1)^— 



1 - 



j.=jV/ ^ '^^nm) Jo 

On using resurgence, i.e., Eq. ([8]), on Tr"*\ namely 
and 



kF 



OO 

r-Z-1 



(r — / — 1)! = / dv exp(— i;) v'^ 
Jo 

after straightforward algebra the second term in Eq. (1A2I) becomes 



/ 1 \ 

u- EE (-1)--^-^ E 



OO ^ 

E TTp 



OO 

X 



X / / — — exp(-M - v) 



Jo u V I 
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or, by interchanging the integral with the r-sum, 

2 oo 



m€An heArn 1=0 {kFnm) 



X / / exp(— M — V) 



Jo u V i_ 



E 



X 

r=AI 



kF 



(A6) 



me,4„ hGAm 1=0 {kFnm) 



m/i 



X / / exp(-M -V, 

Jo Jo U V I '±_ ^ _ / ^ Fnm 

kFnm \u Fmh, 

Eventually, on making the substitution —vF^m/ Fmh ~^ after some algebra the quantity 
in Eq. ( 1A6P can be written as 



_i \N+M 

^ \^ \^ ^7nm+7m?l 

(27ri)2 ^ ^ 



,_„ V ^nm/ 



(AT) 



oo 

X 

1=0 



where the 2nd- level hyperterminant Kn}n{(3; "f) , defined by 



foo POO yTn—l^n—l 

X / / dw df exp(— M — 7^;) ^ r , 

Jo . u\ 



(A8) 
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has been introduced. Accordingly, the complete expression of the 2nd-level hyperasymptotic 
expansion is 

7V-1 
r=0 



1 W 

^ y (-1) 



m£A„ 



M-1 



X 



r=0 



(A9) 



N+M 



(27ri)2 



E E (-1: 



7nm+7 mfi 



n 



r=0 



r,N-M \ kFnm] 



mh 



which coincides with Eq. (1141) . 



APPENDIX B: DERIVATION OF THE EXTRA TERM FOR /? < 



We first recast the double integral in Eq. (fT5|) as follows: 

1 /""^ 



n-l 



Jo ^{u + p){u + v)' 

which, if (3 were not negative, would take the form 

Ki%{(3; 7) = dv exp(-7t;) v 



X 



exp(t;)E^+i(f) - exp(/3)-E^+i(/j) 
(3-v 



(Bl) 



(B2) 
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where the integrand turns out to be continuous even for real, positive values of (3, since 

exp{v)Em+i{v) - exp{P)Em+i{P) 



lim 



(3-v 

1 - exp(/3);3™(m + (3)T{-m, (3) 



(B3) 



If P < 0, Eq. (IBip must be written as 

1 /"°° 



(B4) 



u 

xV I du exp(-u) -, 



where V ... denotes the principal part operator. On evaluating the integral in Eq. (1B4I1 



we have 



u 

V I du exp(-u) -rr-. r 



ml 



exp{v)E,, 



m+ll 



exp(/3)E^+i(/5) 



(B5) 



-ivr exp(/3) 



-py 



p-v 



Finally, the direct substitution of Eq. fIBSP into Eq. flB4p leads to the closed-form of the 
extra term as 



inexp{p){-pr r . , , v^-' 

/ dv expf— 7t>) 

Jo V- p 



P 



n+m—l 



(B6) 
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APPENDIX C: ASYMPTOTIC COEFFICIENTS Tr' FOR THE SWALLOWTAIL 
FUNCTION 

The swallowtail function in Eq. ([T]) corresponds to let g = 1, k = 1, and f{s) = 
— i(s^/5 + xs^/3 + ys'^/2 + zs) in Eq. (1371) . The (four) saddle points are solutions of 
the algebraic equation 

s^ + xs'^ + ys + z = 0, (CI) 
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which can be exactly solved by using, for instance, Cardano's formula. The evaluation of 
the expanding coefficients Tr"'^ requires to evaluate the integral in Eq. ([6]). The first step is 
to expand /(s) — /„ around i.e., 

f{s) -fn = 



- 



+ X ■ 



- si 



+ zis — Sr 



(C2) 



where 



01 



P(s) = + 2SnS^ + 34 + -X S + Us 



Substitution into Eq. ^ gives 



T 



(n) 



WSr 



-X + -y 



(C3) 



1 



(r-1/2)! (5i) 
27rr 



ds 



(C4) 



(r-1/2)! (5i)'+i/2 r 



(2r)! tds2- [P(s)]'-+i/2 

Equation ( 104^ can be further simplified by letting m = s — s„ and then by changing u into 
t = u/ (lOs^ + 5s„x + 5?//2)^/^, which yields 

(n)_ (5ir+^/' (r-1/2)! 



X 



(10s3 + 5s^x + 5i//2)5^/3+i/2 



(2r)! 



]2r 



1 



t=0 



where 



a 



{Wsl + 5snX + 5y/2y/^' 
104 + 5X/3 

^ (10s3 +5SnX + 5?//2)2/3' 

Note that Eq. (ICS!) can also be given the alternative following form: 

(5i)r+l/2 1/2)! 



rpin) 



B. 



(r'+l/2), 



(10s3 + bSnX + 52//2)5'-/3+V2 2r 



a,/3). 



(C5) 



(C6) 



(C7) 
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where the function B^^\u,v) is defined through 

oo ^ 

y rBi^\u,v) = — -— -, (C8) 



n=0 



which coincides with Eq. ( I40l) . 



APPENDIX D: DERIVATION OF THE RECURRENCE RULE IN EQ. (HT]> 

The starting point is the definition of B^\u,v) through the generating function in Eq. 
(]C8|) which, once derived with respect to t, gives 

Et-nBy'K..) = -A (Dl) 

n=0 



and, by taking Eq. ( ICSP into account once again, leads to 

oo 



n=0 



(D2) 



+ (t3 + Mt2^?;t+ 1) ^ = 0. 

n=0 

By operating the products term by term, by rearranging the series indices, after long but 
straightforward algebra Eq. (HTl) follows. 
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LIST OF TABLES 



WT order 


Contribution of saddle S2 


2 


0.02315166515+0.04009986036 i 


3 


0.02303913417+0.03990495099 i 


4 


0.02305968041+0.03994053805 i 


5 


0.02305839073+0.03993830429 i 


6 


0.02305829942+0.03993814613 i 


7 


0.02305830901+0.03993816279 i 


8 


0.02305831068+0.03993816566 i 


9 


0.02305831068+0.03993816562 i 


10 


0.02305831064+0.03993816557 i 


11 


0.02305831064+0.03993816557 i 



TABLE L Values, provided by the WT, of the contribution coming from the saddle S2 in the 
evaluation of the Airy function across its Stokes set, when F = 2. 
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M 


H-WT estimate 


3 


07906105793-0.4032083434 i 


4 


0.7061079570-0.3544207316 i 


5 


0.7018070334-0.3519375922 i 


6 


0.7015955837-0.3518155117 i 


7 


0.7015834126-0.3518084847 i 


8 


0.7015826207-0.3518080275 i 


9 


0.7015825835-0.3518080060 i 


10 


0.7015825939-0.3518080120 i 


11 


0.7015826051-0.3518080185 i 


12 


0.7015826244-0.3518080296 i 


13 


0.7015563466-0.3517928581 i 


14 


0.7015826472-0.3518080428 i 


exact 


0.7015826047 - 0.3518080182 i 



TABLE II: Estimates, provided by the 2nd-level H-WT, of the Airy function across its Stokes set, 
when F = 2. Note that the truncation parameter has been fixed to 15 (corresponding, from 
Fig. [SI to the optimal setting). 
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LIST OF FIGURE CAPTIONS 



FIG. 1: Behavior of the relative error, obtained through the Ist-level H-WT (dots), versus the 
values of A^, which are reported on the abscissa axis. For each value of A^, the values of the relative 
error obtained via the 2nd-level H-WT, with M G [3, — 1], are also plotted and, for the sake of 
clarity, are joined with lines of different color, each of them corresponding to a different value of 
N, departing from the abscissa N itself. 

FIG. 2: The same as in Fig. H but for F = 14 (a), 10 (b), 6 (c), and 2 (d). 

FIG. 3: Behaviors, as functions of F € [2,4], of the relative errors for the Airy function obtained 
via the 1st- (open circles) and the 2nd-level (dots) H-WT. Each point has been obtained by carrying 
out, for a given F, an exhaustive search for those values of the truncations N and {N,M) which 
minimize the corresponding 1st- and 2nd-level relative errors, respectively. 

FIG. 4: Behavior, as a function of F, of the values of corresponding to the optimal setting for 
the Ist-level H-WT in the experiment of Fig. [3l 

FIG. 5: Behavior, as a function of F, of the values of A^ (a) and M (b) corresponding to the 
optimal setting for the 2nd-level H-WT in the experiment of Fig. [3l 

FIG. 6: The same as in Fig. [H but for the instanton integral AA(l/2). 

FIG. 7: Behaviors, as functions of /c G [1/2, 3], of the relative errors for the instanton integral M{k) 
obtained via the 1st- (open circles) and the 2nd-level (dots) H-WT. As for the results presented 
in Fig. [31 each point has been obtained by carrying out, for a given k, an exhaustive search for 
those values of the truncations A^ and (A^, M) which minimize the corresponding 1st- and 2nd-level 
relative errors, respectively. 
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FIG. 8: Behavior, as a function of A;, of the values of N corresponding to the optimal setting for 
the Ist-level H-WT in the experiment of Fig. [71 

FIG. 9: Behavior, as a function of k, of the values of N (a) and M (b) corresponding to the 
optimal setting for the 2nd-level H-WT in the experiment of Fig. [71 

FIG. 10: Behavior of the relative error obtained for the Airy function (dots) and for the instanton 
function (solid curve) versus N, for F = k = 3 (a), 7 (b), 12 (c), and 20 (d). 

FIG. 11: Pictorial representation of the saddle networks and of the complex integration paths 
involved in the evaluation of the Airy (a) and of the instanton (b) functions. 

FIG. 12: The same as in Fig. [TTl but for the evaluation of the swallowtail diffraction catastrophe 
at the points across the Stokes set given by the triplets (x, y, z) = (0, k^/^, x 0.23012 . . .), with 

K>0. 

FIG. 13: The same as in Fig. [H but for the swallowtail function evaluated across the Stokes set 
(x, y, z) = (0, k2 X 0.23012 . . .) with k = 2. 

FIG. 14: Behaviors, as functions of k G [2,4], of the relative errors for the swallowtail function 
evaluated across the Stokes set {x,y,z) = (0,k;^/^,k^ x 0.23012...) obtained via the 1st- (open 
circles) and the 2nd-level (dots) H-WT. As for the results presented in Fig. [3l each point has been 
obtained by carrying out, for a given k, an exhaustive search for those values of the truncations 
and (A^, M) which minimize the corresponding 1st- and 2nd-level relative errors, respectively. 

FIG. 15: Behavior, as a function of k, of the values of A^ corresponding to the optimal setting for 
the Ist-level H-WT in the experiment of Fig. [TH 

FIG. 16: Behavior, as a function of k, of the values of A^ (a) and M (b) corresponding to the 
optimal setting for the 2nd-level H-WT in the experiment of Fig. [TH 



33 




34 




35 




36 



Optimal N 
10 



8 



O 



CHD GQQOQOQQOQQOOQOQO 
QOOOOOOOOOOOOOOOOOOOOOOO GOO 

4mMm (Mm 



, , , , \ , , , , \ , , , , \ , , , , L p 

2.0 2.5 3.0 3.5 4.0 

Fig. 4 - R. Borghi 



37 



(a) 

Optimal N 
30 



25 
20 
15( 
10 



cr" OO ^ 
(B^ O® O OO 



2.0 2.5 3.0 3.5 4.0 

(b) 

Optimal M 
30 

25 

20 

15 



10^ 



5 



^QX ^ O 



O CSX o - 

C59 O 



2.0 2.5 3.0 3.5 4.0 

Fig. 5 - R. Borghi 



38 



Relative error 




Fig. 6 - R. Borghi 



39 



Relative error 




Fig. 7 - R. Borghi 



40 



Optimal N 

10 r 



8 - 

6 - ceo 

000000 0000000000 
4<300O 0000000000000000000000000 000 

2 - 

0.5 1.0 1.5 2.0 2.5 3.0 

Fig. 8 - R. Borghi 



41 



(a) 



Optimal N 
30 



25 
20 
15 
10 

5' 



o ° ^ 

O ^00 



O GOD qO 

^ GO 



0.5 1.0 1.5 2.0 2.5 3.0 

(b) 

Optimal M 
30 

25 

20 

15 

10 

5' 



O O O 



o 



0.5 1.0 1.5 2.0 2.5 3.0 

Fig. 9 - R. Borghi 



42 




43 



(a) 




O Si 



Fig. 12 - R. Borghi 



45 



Relative error 




Fig. 13 - R. Borghi 



46 



Relative error 




2.5 3.0 3.5 4.0 

Fig. 14 - R. Borghi 



47 



Optimal N 
30 



25 
20 
15 
10 

5. 



OOO^ 
OOOO^ 

oooo^ 
oooo^ 

00000°°° 



— I 1 1 1 \ 1 1 1 1 \ 1 1 1 1 \ 1 1 1 1 L ^ 

2.0 2.5 3.0 3.5 4.0 

Fig. 15 - R. Borghi 



48 



(a) 

Optimal N 
30 

25 

20 

00000000000000000000 

15 ; 

0000000000000000000 

5CED 



^ , , , 1 , , , , 1 , , , , 1 , , , , 1_ 

2.0 2.5 3.0 3.5 4.0 



(b) 

Optimal M 
30 

25 

20 

15 

GOOOOOOOOOOOOOOOOOOO 

10 

' 0000000000000000000 

; ; ; 1 ; ; ; , 1 ; , ; ; 1 ; ; , ; L 

2.0 2.5 3.0 3.5 4.0 

Fig. 16 - R. Borghi 



49 



